On the two-proton emission of Fe - a new type of radioactivity 
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In an experiment at the SISSI-LISE3 facility of GANIL, 
the decay of the proton drip-line nucleus 45 Fe has been studied 
after projectile fragmentation of a 58 Ni primary beam at 75 
MeV/nucleon impinging on a natural nickel target. Fragment- 
implantation events have been correlated with radioactive de- 
cay events in a 16x16 pixel silicon strip detector on an event- 
by-event basis. The decay-energy spectrum of 45 Fe implants 
shows a distinct peak at (1.06±0.04) MeV with a half-life of 
T1/2 = (4.7^j' 4 ) ms. None of the events in this peak is in coin- 
cidence with f3 particles which were searched for in a detector 
next to the implantation detector. For a longer correlation 
interval, daughter decays of the two-proton daughter 43 Cr 
can be observed after 45 Fe implantation. The decay energy 
for 45 Fe agrees nicely with several theoretical predictions for 
two-proton emission. Barrier-penetration calculations slightly 
favour a di-proton emission picture over an emission of two 
individual protons and point thus to a 2 He emission mode. 

PACS numbers: 23.50.+Z, 23.90. +w, 21.10.-k, 27.40. +z 

An ensemble of protons and neutrons can form a nu- 
cleus stable against any radioactive decay only if a subtle 
equilibrium between the number of protons and neutrons 
is respected. If this equilibrium condition is not respected 
in a nucleus, it becomes radioactive. For small deviations 
from equilibrium, the nuclei decay by [3 decay. The lim- 
its of stability, the drip lines, are reached if the nuclear 
forces are no longer able to bind an ensemble of nucleons 
with a too large neutron or proton excess. Whereas for 
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heavy nuclei, e.g. a decay, 14 C radioactivity or fission 
occurs, unstable proton-rich nuclei may decay by emis- 
sion of one proton for odd-Z nuclei or of two protons for 
even-Z nuclei from their ground states. The one-proton 
ground-state decay has been observed for the first time 
at GSI Darmstadt in 1981 @J|. Meanwhile almost 30 
cases of proton radioactivity have been identified for odd- 
Z nuclei beyond the proton drip line between Z=51 and 
Z=83 H allowing e.g. to establish the sequence of shell- 
model single-particle levels beyond the proton drip line. 

Two-proton (2p) radioactivity is predicted since 
1960 B to occur for even-Z proton-rich nuclei beyond the 
drip line. Due to the pairing energy, the 2p candidates 
cannot decay by a sequential emission of two protons as 
the one-proton daughter is energetically not accessible. 
Therefore, only a simultaneous two-proton emission is 
possible which can take place in two different ways: i) 
by an isotropic emission of the two protons which then 
have no angular correlation, i.e. they fill the whole phase 
space available, but in order to easily penetrate through 
the Coulomb and centrifugal barrier of the daughter nu- 
cleus they share most probably equally the decay energy; 
ii) by a correlated emission where in the decay the 2 He 
resonance is formed which decays either already under 
the Coulomb and centrifugal barrier or outside the bar- 
rier. In both cases, a zero energy difference between the 
two protons is most likely. However, for a 2 He emission, 
a small relative angle between the two protons might be 
observable. 

Recent theoretical work [|J-|7) showed that 45 Fe, 48 Ni, 
and 54 Zn are the best candidates for two-proton ground- 
state decay as their 2p Q values are about 1.1-1.8 MeV, 
whereas one-proton emission is either energetically for- 
bidden or extremely disfavoured due to small one-proton 
decay energies and very narrow intermediate states. 

For light nuclei, the Coulomb and centrifugal barri- 
ers for proton emission are smaller and the width of the 
ground states of these nuclei become large. Thus, onc- 
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proton emission is energetically possible through the tails 
of broad intermediate states even for even-Z nuclei. This 
has been observed experimentally for 6 Be [|| and 12 O 0. 
In both cases, the emission pattern is compatible with a 
simultaneous uncorrelated decay into the available phase 
space. Two-proton emission from excited states has been 
observed in eight cases after (3 decay (see JlO| for a recent 
reference) as well as after inelastic excitation of 14 O |Tl| l 
and of 18 Ne 12 1. In all cases, the experimental data are 
either not precise enough to indicate a detailed decay 
pattern or compatible with a sequential or three-body 
picture with half-lives as short as 10~ 21 s. 

With the advent of projectile-fragmentation facilities 
equipped with powerful separators for in-flight isotope 
separation, medium-mass proton drip-line nuclei came 
into experimental reach and two of the above mentioned 
most promising candidates, 45 Fe Jl3| and 48 Ni [[uj, could 
be observed. In addition, part of the decay strength 
of 45 Fe was observed 15 in the experiment designed to 
identify for the first time 48 Ni iQ. However, due to the 
set-up of the electronics in this experiment, two-proton 
events triggered the data acquisition only with a very low 
probability (five possible events at about 1.1 MeV 
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FIG. 1. Two-dimensional fragment identification spectrum 
for the present experiment. The implantation events are plot- 
ted as a function of their time-of-flight between one chan- 
nel-plate detector and the silicon stack and their energy loss 
in the first silicon detector. The figure shows only part of the 
data. 

In an experiment performed in June/July 2000 
at the SISSI-LISE3 facility of GANIL, we used the 
projectile fragmentation of a 58 Ni primary beam at 
75 MeV/nucleon to produce proton- rich nuclei in the 
range Z=20-28. After production in a Tmt Ni target (240 
/im in thickness) located in the SISSI device, the frag- 
ments of interest were selected by the Alpha/LISE3 sep- 
arator equipped with an intermediate beryllium degrader 
(50/im). At the focus of the LISE3 separator, a set- 
up was mounted to identify (see figure |l|) and stop the 
fragments as well as to study their radioactive decays. 



This set-up consisted in two channel-plate detection sys- 
tems for timing purposes mounted at a first LISE focal 
point 22.9 m upstream from the final focus, a sequence 
of four silicon detectors (El: 300^m, E2: 300/im, E3: 
300/im, E4: 6mm, respectively) with the third one being 
a silicon-strip detector with 16x16 x-y strips with a pitch 
of 3 mm, and a Germanium array in close geometry. The 
silicon detectors were equipped with two parallel elec- 
tronic chains with different gains, one for heavy- fragment 
identification and the other for decay spectroscopy. 

The fragments of interest were stopped in the third sil- 
icon detector of the telescope and identified on an event- 
by-event basis by means of their flight times and of their 
energy loss in all detectors of the telescope (see Jl4|,|l5| ) ■ 
Implantation events were triggered by detectors El and 
E2, whereas radioactive decay events were triggered ei- 
ther by the detector E3 or by the adjacent silicon detector 
E4. The efficiency to observe a (3 particle in the E4 detec- 
tor for a (3 decay occurring in the implantation detector 
is about 30%. The 7 energy calibration and detection ef- 
ficiency was obtained with standard calibration sources. 
The efficiency was about 1.6% at 1.3 MeV. 

In a 36 h run, 22 45 Fe implantations were identi- 
fied. Due to a rather low implantation rate of much less 
than one radioactive isotope per second in each pixel, an 
implantation-decay correlation could be performed on an 
event-by-event basis with each decay event being corre- 
lated to each implantation in the same x-y pixel having 
occurred less than a well-defined maximum time interval 
before the decay event. 
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FIG. 2. Decay-energy spectrum correlated with 45 Fe im- 
plantation. Spectrum a) is obtained under the condition that 
the radioactive decay occurs faster than 15 ms after implan- 
tation, whereas spectrum b) contains decay events with decay 
times between 15 ms and 100 ms. The peak at 1.06 MeV is 
clearly identified as being due to a fast decay. 

Figure ^a shows the decay-energy spectrum correlated 
with implants of 45 Fe where only decay events occur- 
ring less that 15 ms after a 45 Fe implantation were anal- 
ysed. The spectrum exhibits a pronounced peak at 
(1.06±0.04) MeV with only a very few other counts. In 
contrary, in the spectrum in figure |2|b conditioned by a 
decay time in the interval between 15 ms and 100 ms, 
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the 1.06 MeV peak has almost completely disappeared 
and other events higher in decay energy show up. These 
counts are consistent with the decay-energy spectrum of 
43 Cr ||, the 2p daughter of 45 Fe. The 1.06 MeV peak, 
however, seems to originate only from the fast decay of 
45 Fe. In addition, the events in this peak have no coin- 
cident /3-particle signals in the adjacent detector E4 (see 
figure |3|a) beyond the noise level, whereas these coinci- 
dent (3 particles can be observed for the (3 decay of 46 Fe 
implants analysed with a similar condition in energy for 
the E3 detector. For 12 46 Fe events in the decay-energy 
interval 0.875 - 1.3 MeV, six coincident (3 particles can 
be identified in the E4 detector (see figure ||b). As the 
/3-decay end-point energies are roughly the same for 45 Fe 
and 46 Fe and therefore similar /3-detection efficiencies can 
be assumed, the non-observation of (3 particles in coinci- 
dence with the 1.06 MeV peak alone is a strong indica- 
tion that this peak originates from a direct two-proton 
ground-state decay of 45 Fe. The probability to miss all 
(3 particles for the 12 events in the peak is as low as 
1.4%. Although the 7-detection efficiency was rather low, 
it is worth mentioning that none of the 12 events in the 
1.06 MeV peak is in coincidence with a 7 ray. 
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FIG. 3. a) /3-particle spectrum from the 6mm-thick detec- 
tor E4 in coincidence with events in the peak at 1.06 MeV. 
b) Similar spectrum obtained after 46 Fe implantation con- 
ditioned by a decay-energy range in the E3 detector of 
E — 0.875 — 1.3 MeV, which yields the same statistics as for 
the 1.06 MeV peak of 45 Fe. The insets show the decay-energy 
spectrum for the two nuclei with the dark area and the dotted 
lines indicating the decay-energy condition applied to gener- 
ate the (3 spectra as well as the full-statistics 46 Fe (3 spectrum. 

The observation of only twelve events in the peak at 
1.06 MeV indicates that the branching ratio for 2p decay 
of 45 Fe is not hundred percent, in agreement with theo- 
retical prediction of a /3-decay half-life of about 7 ms || . 
However, dead-time losses (the data-acquisition dead 
time is about 0.5 ms) make us lose about 3-4 events. 
The dead zone between two strips of the implantation 
detector is another possible source of losses. The events 
with an energy above 6 MeV (see figure ||a) decay with 
a half-life compatible with the one of 45 Fe and are there- 
fore must likely /3-delayed events. With this information, 
a branching ratio for 2p emission of 70-80% can be esti- 



mated. 

The observation of a pronounced peak at 1.06 MeV 
from the decay of 45 Fe is a strong additional indica- 
tion that 45 Fe is not decaying by a /3-delayed mode 
alone as does 46 Fe |l5| with a more complexe decay- 
energy spectrum. In addition, the width of the Fe 
peak (60±10 keV) is about 40% smaller than /3-delayed 
one-proton peaks from neighboring nuclei. Although the 
statsistics is limited, this indicates that the 45 Fe peak is 
not broadened due to /3 pile-up. 

The decay-time spectrum of 45 Fe gated by the 
1.06 MeV peak is shown in figure || A one-component fit 
with an exponential yields a half- life for 45 Fe of TV2 — 
(4.7i.f; 4 ) ms. The decay-time spectrum of events up to 
100 ms after a 45 Fe implantation can be fitted by taking 
into account the decay of 45 Fe and its 2p daughter 43 Cr. 
The half-life then is (5.7±?;^) ms. 

The energy of the peak at 1.06 MeV agrees nicely with 
Q2 P value predictions from Brown || of 1.15(9) MeV, 
from Ormand § of 1.28(18) MeV, and from Cole § 
of 1.22(5) MeV. These models use the isobaric-multiplet 
mass equation and shell-model calculations to determine 
masses of proton-rich nuclei from their neutron-rich mir- 
ror partners. Q-value predictions from models aiming 
at predicting masses for the entire chart of nuclei are in 
less good agreement with our present result. All the 2p 
Q-value predictions are summarized in figure ||, where 
they are used in barrier-penetration calculations using 
the simple di-proton model for two-proton emission (the 
model used in H with R = 4.2 fm). 
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FIG. 4. Decay-time spectrum of 45 Fe. The twelve events in 
the 1.06 MeV peak are fitted by a one-component exponential 
using the maximum likelihood procedure. The half-life thus 
obtained is (4.7t.i't) ms - A fit, including the daughter decay, 
of all correlated event up to 100 ms after 45 Fe implantation 
(see inset) yields a half-life of (5.7if] 4 ) ms for 45 Fe and a 
value consistent with literature fla] for 43 Cr. The dark shaded 
counts originate from the 1.06 MeV peak. 

For E p = 1.06±0.04 MeV, these calculations predict a 
di-proton barrier penetration half-life of (0.21q 3 ) ms, if 
one assumes a spectroscopic factor of unity. Brown pjlq] 
calculated a spectroscopic factor of 0.195 for a direct 2p 
decay of 45 Fe which increases the barrier tunnel time to 
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a predicted half-life of (1.2 jl g) ms. 

We have also used the R-matrix formulation of 
Barker p3| both in its simple form, equivalent to the 
di-proton decay model discussed above, and in a more 
realistic form which includes the s-wave resonance of 
the two protons. In the R-matrix we need to choose 
a channel radius a and an associated dimensionless re- 
duced width # 2 p (Eq. 16 of For consistency we 
take a = 4.2 fm which is a region of the resonance wave 
function where 2 w 1. With S = 0.195 we obtain a 

half-life of (1.7to'|) ms m tne simple model, and (8t;P) s 
in the s-wave resonance model. The reduction in the 
phase space due to the two-proton resonance increases 
the lifetime by a factor of 4000. Grigorenko et al. ]25j] 
studied the emission of two protons from 48 Ni in a three- 
body model. They find that the half-life for the same 
Q 2p value increases by about three orders of magnitude 
when going from a simple di-proton model to the more 
advanced three-body model. Thus the models which in- 
clude the interaction between two-protons appear to give 
a lifetime which is much longer than the observed value. 

A sequential decay seems to be excluded for 45 Fe, as 
the intermediate state, the ground state of 44 Mn, is, de- 
pending on the prediction used, either not in the allowed 
region or at the very limit of the allowed region. The 
model predictions (^-0] range from Q\ p = -24 keV to 
+10 keV. As the intermediate state is most probably 
rather narrow (barrier-penetration calculations yield a 
value of about Y = 50 meV), the first proton would have 
a very long half-life (hours or even days), which makes 
this decay mode very unlikely. 
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FIG. 5. Barrier-penetration half-life as a function of the 
two-proton Q value, Qi v , for 45 Fe. The barrier penetra- 
tion was calculated by assuming a spectroscopic factor of 
unity. Different model predictions |5|-|7|,[17H22| were used for 
Qi v . The experimentally observed Q value of 45 Fe implies a 
di-proton barrier-penetration half-life of 0.2 ms. 

In conclusion, we have studied the decay of 45 Fe in an 
experiment where we measured its decay energy yielding 
a peak at 1.06 MeV, a half-life of 4.7 ms, and, in particu- 
lar, no coincident (3 particle or 7 ray for the events in the 
decay-energy peak. In addition, strong indications for 
the decay of the 2p daughter, 43 Cr, after implantation 



of 45 Fe are found. Additional support comes from the 
width of the 45 Fe peak. The energy of the observed peak 
is in reasonable agreement with theoretical predictions 
for the 2p decay energy. A consistent picture arises, if 
one assumes that a two-proton ground-state emission oc- 
curs. Simple model calculations favor the observed decay 
mode to be a 2 He/di-proton emission. 

Whereas the two-proton ground-state decay of 45 Fe 
seems to be established with the present data, future 
high-statistics data should definitively allow to conclude 
on the nature of the two-proton decay, 2 He emission or 
three-body decay. This question can be addressed in an 
experiment, which measures the individual proton ener- 
gies and the relative- angle distribution for the two pro- 
tons emitted which should be either isotropic (three-body 
decay) or forward-peaked ( 2 He emission). 

It should be mentioned that similar results, although 
with less statistics and a lower energy resolution, were 
also obtained at the FRS of GSI §|]. 
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